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The onset and characteristics of Micro- Tearing Modes (MTM) in the core of spher- 
ical (NSTX) and conventional tokamaks (ASDEX-UG and JET) are studied through 
local linear gyrokinetic simulations with gyro [J. Candy and E. Belli, General Atom- 
ics Report GA-A26818 (2011)]. For experimentally relevant core plasma parameters 
in the NSTX and ASDEX-UG tokamaks, in agreement with previous works, we 
find MTMs as the dominant linear instability. Also, for JET-like core parameters 
considered in our study an MTM is found as the most unstable mode. In all these 
plasmas, finite collisionality is needed for MTMs to become unstable and the electron 
temperature gradient is found to be the fundamental drive. However, a significant 
difference is observed in the dependence of linear growth rate of MTMs on electron 
temperature gradient. While it varies weakly and non-monotonically in JET and 
ASDEX-UG plasmas, in NSTX it increases with the electron temperature gradient. 

I. INTRODUCTION 

In recent years, specially in view of an increasing interest in high /3 e operation scenarios, 
such as hybrid scenarios for ITER , the impact of electromagnetic effects on the particle 
and heat transport has attracted much attention. Here /3 e = (8irn e T e ) / B 2 nit , where n e and 
TL are the electron density and temperature. B unit is defined as the effective field strength 
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Recent reports have shown the significant role of electromagnetic modes such as Micro- 
Tearing Modes (MTMs) on the electron heat transport in the core of fusion plasmas 4rtZ|. 
It has been found that in plasmas where /3 e and collisionality are sufficiently high MTMs 
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can become the dominant instability. Under these conditions, MTMs generate a major 
contribution to the anomalous electron heat flux. 

The MTMs are small-scale in the radial direction, but ion-scale in the binormal direction. 
They are electromagnetic modes with an even parity with respect to the perturbed parallel 
vector potential, 8 An (referred to as tearing parity). In the literature two drive mechanisms 
for MTMs are proposed: one is the time-dependent thermal force experienced by electrons 
which results in a parallel current that produces magnetic field perturbations. If these 
perturbations then tip the field lines in the direction of an equilibrium electron temperature 
gradient the thermal force will increase and therefore the instability arises The second 
mechanism is due to a current carried by nearly trapped electrons in a boundary layer close 
to the trapped-passing boundary j^, 1(|. Both of these mechanisms require finite electron 
temperature gradient and collisions. However, MTMs have been observed to arise under 
conditions which are not well-described by the above mechanisms \&, Thus, a complete 
picture of MTM excitation is not available at present. 



In the following we summarize the previous findings reported in Refs. 4|, lllMl4j . Linear 



simulations for MAST 



NSTX[12[, and ASDEX-UG [14J plasmas reported that MTMs 



can be the dominant instabilities in the region r/a = 0.5 — 0.8 with the maximum growth 
rate at mode numbers kgp s between 0.2 — 0.8. Here, p s is the ion sound Larmor radius. The 
electron temperature gradient is found to be the drive of the instability and a non-monotonic 
dependence of the growth rate on the electron-ion collisionality v e i is observed. The peak 
of the growth rate coincides with the experimental value of v e i for various considered radial 
positions. The non-monotonic dependence of the MTM growth rate on collisionality is due 
to the fact that, on one hand both of the driving mechanisms mentioned above require finite 
collisionality, and therefore stability can be expected as v ei reduces, but on the other hand, 
in a strongly collisional regime the strong rate of scattering of the electrons between the 
field lines prevents the formation of a current layer, hence MTMs are stabilized fl^ |. 

It has also been shown that by increasing the effective ion charge Z e // ; MTMs are desta- 
bilized through the Z e j/-dependence of the electron-ion collision frequency, see Ref. jl2 |. 
As MTMs are electromagnetic in nature a finite /3 e is needed for their destabilization. For 
the studied discharge in Ref. 12| it is found that the /3 e threshold is well below the ex- 
perimentally relevant /3 e value and the growth rate increases moderately with increasing 
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Non-linear simulations confirmed the role of the electron temperature gradient as the 



drive of the MTM instability in spherical (NSTX) [5[ as well as standard tokamaks (ASDEX- 
UG) 14j. Moreover, it has been shown that the Chirikov criterion 15[ for overlapping of 
ihe magnetic islands leads to an up-shift of the electron temperature gradient threshold 
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In the present paper we investigate the onset of the MTMs and its parametric dependence 



through local linear gyrokinetic simulations with the GYRO code 



16|, in a spherical tokamak: 



NSTX, and two conventional tokamaks: ASDEX-UG and JET. For the NSTX case we use 



plasma parameters found in Ref. 



the plasma parameters reported in Ref. 12|, LL3[], and for the ASDEX-UG case we use the 



In the present work we re-examine these discharges 



with an emphasis on the parametric dependences of the MTM onset. In view of the new 
ITER-like wall experiments on JET, in the presented analysis we have also considered a set 
of JET-like parameters. 

The remainder of the paper is organized as follows. In Sec. [IT] the input parameters 
are discussed, and in Sec. DTD parametric dependences of the MTM onset are analyzed by 
presenting scans over MTM driving parameters such as collisionality, /3 e and electron tem- 
perature/density scaling lengths. The conclusions are drawn in Sec. IIV1 



II. INPUT PARAMETERS 

The plasma parameters used in our analysis are shown in table [H 



TABLE I: Input parameters for densities, temperatures and their gradients. 





r/a 


Z eff 


n e [10 19 /m 3 } 


T e [keV] 


a /L ne 


a/L Ti 


a/L Te 


Ti/T e 


Ps/a 


v ei {a/c s ) 


NSTX 


0.6 


2.91 


6.0 


0.44 


-0.83 


2.36 


2.72 


0.94 


0.0074 


1.45 


AUG 


0.65 


3.30 


7.6 


0.765 


0.37 


2.18 


3.02 


1.19 


0.0018 


0.68 


JET 


0.6 


3.41 


7.8 


1.25 


0.15 


2.16 


2.16 


1.00 


0.0027 


0.43 



Here, L n = — [<9(lnn)/<9r] _1 , Lj< = — [<9(lnT)/<9r] -1 , are the density and temperature 
scale lengths, a is the outermost minor radius. Deuterium ions, an active impurity species 
denoted in table [IT] by Z\ (carbon for NSTX and nitrogen for ASDEX-UG and JET), and a 
passive species of impurity (tungsten unless otherwise stated) denoted in the table [TT1 by Z 2 , 
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TABLE II: Input parameters for plasma ion compositions. 







nzi/n e 


a/L n zi 


z 2 


nz2/n e 


NSTX 


C +6 


6.4% 


-2.75 


w +m 


0.02% 


AUG 


N+ 7 


4.8% 


0.80 


w +m 


0.02% 


JET 


N+ 7 


5.0% 


0.14 


w+ 40 


0.02% 



are considered. The passive species are considered here in order to examine the impurity 
particle transport due to MTMs. Note that in the ASDEX-UG case reported in Ref. 14] 
no impurities were present and the value of a/Lxi has been artificially reduced from its 
experimental value to eliminate the drive of Ion Temperature Gradient (ITG) modes, but 
here we use the experimental values. 

Linear runs with gyro include full electromagnetic effects: shear SB (= V x 5 Ah), 
and compressional 8Bn magnetic perturbations. Gyrokinetic electrons are assumed, and 
the collisions are modeled using an energy dependent Lorentz operator. Both electron-ion 
and electron-electron collisions are included in the electron collision frequency z/ e (w), and 
collisions between all ion species are accounted for. To take into account the plasma shape we 
have used a Miller- type local equilibrium model available in gyro, see Refs. 3,[l7|. Typical 
JET parameters for plasma shape and magnetic geometry are used, and the corresponding 
values are given in Table IIIIl In this table, (3 e is calculated in CGS units following the 



TABLE III: Input parameters for plasma shape and magnetic geometry. 





a[m] 




OiMHD 


Bunit 


R /a 


Bo[T] 


s 


q 


K 


<5 


kep s 


NSTX 


0.6 


0.024 


0.36 


0.66 


1.52 


0.35 


1.73 


1.68 


1.72 


0.12 


0.63 


AUG 


0.6 


0.005 


0.42 


2.16 


3.3 


2.479 


1.31 


2.18 


1.30 


0.13 


0.2 


JET 


1.0 


0.013 


0.35 


1.77 


3.3 


2.55 


1.32 


1.45 


1.70 


0.37 


0.5 



expression: 



8vr(n e [10 19 /m 3 ] 1(T 6 10 19 )(T e [keV] 1.6022 10" 



(1) 



(10 4 B wU [T\Y 

where B un u = (d'Xt/dr )/r is the effective field strength with \t defined through the toroidal 
magnetic flux 27r\; t jsllol, q is the safety factor, and the magnetic shear is s = (r/q)dq/dr. 
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The generalized magnetohydrodynamic a parameter is defined as 

2d 8yr d P (o\ 
0(mhd - ~q K —^——c p , {2) 

D unit ar 

where Rq is the effective major radius, r is the minor radius, and p = Y2 a n aT a is the total 
plasma pressure. c p is a scaling parameter which allows an artificial adjustment of oimhd 
(affecting the magnetic curvature drift) without modifying the background gradients, as 
presented in Ref. 18[]. Furthermore, k is the elongation, 5 is the triangularity. 

Typical resolution parameters used in our linear analysis are as follows: 40 radial grid 
points, 12 parallel orbit mesh points (x2 signs of parallel velocity), 16 pitch angles, and 8 
energies. A high radial resolution is needed since the linear instability of MTM depends 
on the presence of narrow resonant current layers centered on the rational surfaces. The 
toroidal mode numbers used here, corresponding to the p s /a values in table [H are n = 30 
for NSTX, n = 18 for AUG, and n = 75 for JET. 



III. MTM INSTABILITY 

Figures[T](a,b) illustrate the linear growth rate and the real frequency of the most unstable 
modes in the three machines as functions of kgp s . In NSTX and JET, the MTMs (with 
positive real frequency in electron diamagnetic direction) are the most unstable modes in 
the range of 0.1 < kgp s < 1, while for ASDEX-UG the MTMs are the most unstable 
modes only in the narrower wave number range 0.1 < kgp s < 0.3, corresponding to longer 
wavelengths than those in NSTX and JET. In NSTX the MTMs remain the most unstable 
mode for the whole considered range of kgp s , but in ASDEX-UG an ITG mode (negative 
real frequency, in ion diamagnetic direction) for 0.3 < kgp s < 1 and an ETG mode (positive 
real frequency, in electron diamagnetic direction) for 1 < kgp s are also found to be unstable. 
The kgp s corresponding to the maximum growth rate varies between the different machines: 
in NSTX k e p s ~ 0.6, in ASDEX-UG k e p s ~ 0.2 and in the JET like case k e p s ~ 0.5. 
The normalized poloidal mode number for the remainder of our calculations are fixed to 
these values (see table ITTTT) . Previous studies have discussed the difference in the poloidal 
mode number corresponding to the maximum of the unstable MTMs between the NSTX 
and ASDEX-UG cases, and it is believed to be due to characteristics of the spherical or 
conventional tokamaks [12I . l^ . However, here we find unstable MTMs with similar mode 
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FIG. 1: Imaginary (a), and real parts (b) of eigenvalues (7, u r ) as functions of kgp s . Red solid 
lines: NSTX, blue dashed: ASDEX-UG, green dash-dotted: JET. 



numbers in the JET tokamak to that of NSTX, which suggests that this is not always the 
case. Therefore, further studies are needed to determine the reason for the similarities and 
differences in the mode numbers for maximum growth rates in various machines. 

The structure of the 8cj), SA\\ and 8B\\ eigenmodes, corresponding to the kgp s values 
mentioned above are shown in figures [2] (a-f). The eigenfunctions are normalized so that 
8A\\ [9 = 0) is unity. The MTM signature is distinguished by the tearing parity of the 8A\\ 
eigenmodes. The NSTX case has the strongest electromagnetic character in terms of the 
relative amplitude of 8<j), while the "most electrostatic" mode is the AUG case. Regarding 
the strength of the compressional magnetic perturbations, 5B\\ , the situation is the opposite; 
the less electrostatic the mode is the stronger the 8B\\ perturbations are. The characteristic 
width of the eigenmodes along the field line is similar in all three machines. 

Figures |3] and H] illustrate the normalized linear energy and particle fluxes, respectively, 
and the contributions to these from the 8<p, 8A\\ and 8B\\ fluctuations, as functions of kgp s 
in all three machines. As seen in figures [3] (a-c), in the MTM dominated region the electron 
heat flux is the main channel of the energy transport, and the dominant contribution is 
generated by 8A\\ (blue dash-dotted line). In the NSTX case, see figure [3] (a), this is true 
for the whole considered range of k$p s . The ion heat flux generated by MTMs shown in 
figure |3] (d) (magenta solid line) however, is negligibly small compared to the electron heat 
flux. The particle fluxes generated by MTMs for electrons as seen in figure H] (a-c), and both 
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FIG. 2: Linear parallel mode structures of 5<f> (a,b), 8 An (c,d) and SBu (e,f). Note that the actual 
radial resolution of the simulations covers 9/ir = (—16, 16). 



active and passive impurity species shown in figures H] (d-f and g-i) are also negligible in 
comparison to the electron heat flux. We note, that nonlinear MTM simulations for NSTX, 
presented in js[ is| also showed negligibly small particle fluxes. 

In the case of ASDEX-UG the main contribution to electron heat flux, see figure [3] (b), is 
generated by the MTM instability at low kep s , and for higher kop s where the most unstable 
mode switches to an ITG mode, the electron heat flux is significantly reduced. However, at 
the very high poloidal mode numbers 1 < kgp s where an ETG is the dominant instability 
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NSTX 



AUG 




FIG. 3: Normalized linear electron (top) and ion (bottom) heat fluxes (magenta solid lines) and 
their contributions from S(f> (red dashed lines), SA^ (blue dash-dotted lines), and 5B\\ (green dotted 
lines) versus k ePs in NSTX (a, d), ASDEX-UG (b, e), and JET (c, f). 



the electron heat flux increases again. Also here, 5A\\ (blue dash-dotted line) generates the 
dominant contribution to the electron heat flux in the MTM dominated region while 5<p (red 
dashed line) produces the dominant contribution to the electron heat flux in the ITG/ETG 
dominated regions. For the ion heat flux the main contribution comes from the higher kgp s 
region where the ITG mode is the most unstable mode present with the maximum around 
kgp s ~ 0.5, as illustrated in figure [3] (e). The MTM and ETG contributions to the ion heat 
flux are significantly smaller. Also, for the electron and impurity particle fluxes, shown in 
figures H] (b, e and h), the contributions from MTM and ETG is negligible compared to the 
contribution from the ITG. 

Similar trends are observed for heat and particle transport in the JET case. Again here, 
in the MTM dominated region the electron heat flux, see figure [3] (c), is the main channel 
of transport while the ion heat flux, presented in figure [3] (f), and the particle fluxes shown 
in figures H] (c, f and i) are negligibly small in comparison. 

Remarkably, the main ion energy fluxes (Qi) generated by the ITG modes in the ASDEX- 
UG is found to be inward in spite of the positive ion temperature gradient, see Figs. [3] (d). 
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FIG. 4: Normalized linear particle fluxes (magenta solid lines) and their contributions from 6(j> (red 
dashed lines), 8 An (blue dash-dotted lines), and SBi\ (green dotted lines) versus kep s in NSTX (a, 
d and g), ASDEX-UG (b, e and h), and JET (c, f and i). 



However we note, that in the ASDEX-UG case 1) there is a high inward particle flux of 
main ions which might account for the inward energy flux if most of it is convective, 2) the 
strong positive ion energy flux carried by the nitrogen impurities almost cancel the negative 
energy flux of the main ions, so the total ion energy flux is close to zero. 



A. Parametric dependences 

Figures [5] (a,c) show the linear growth rates and the real frequencies of the unstable 
modes, corresponding to the fixed values of k$p s given in table [TTT| as functions of j3 e for the 
three considered machines. As seen in this figure for the experimental values of f3 e (shown 
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with vertical lines) the most unstable mode is an MTM for all three tokamaks. In NSTX, 
the onset of MTM is well below the experimental value of /3 e and an increase in /3 e above its 
experimental value does not increase the growth rate significantly. The same observation can 
be made for the JET case, while for ASDEX-UG, although the MTM onset is well below its 
experimental value, by further increasing /3 e , a Kinetic-Ballooning Mode (KBM) (negative 
real frequency in ion direction) appears as the most unstable mode. These results show that 
since MTMs become unstable due to electromagnetic perturbations a finite level of (3 e is 
needed for their onset, but when unstable, MTMs do not strongly depend on (3 e variations. 
In all these scans the value of olmhd is calculated consistently with the local beta values 
and the density and temperature gradients. 
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FIG. 5: Imaginary and real parts of eigenvalues (7, uj r ) as functions of f3 e (a and c), and v e i (b and 
d). Red solid lines: NSTX, blue dashed: ASDEX-UG, green dash-dotted: JET. The vertical lines 
represent the base parameters color coded similarly to their respective machine. 



In order to determine the effect of olmhd through the curvature drift a /3 e scan, similar to 
figure [S] (a and c), is performed where olmhd is scaled to zero by setting c p = 0. The corre- 
sponding results are shown in figure [6j The «M#D-stabilization is not significant in the MTM 
regime in any of the studied plasmas. In the NSTX case there is a small ctiMfrD-stabilization 
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only at very high j3 e , see figures [6] (a,d). In the ASDEX-UG case, the stabilization of the 
KBM mode at higher /3 e is clearly an omhd effect, since without it the mode is further 
destabilized by an increase in /3 e , see the black dashed lines in figures [6] (b,e). For JET, the 
situation is different, as seen in figures |6] (c,f). In the electrostatic limit and at low /3 e the 
most unstable mode is an ITG which is stabilized as /3 e is increased even without oimhd 
effect. For higher /3 e the dominant mode switches to an MTM. Without oimhd effects, by 
increasing f} e even further the MTM switches to a KBM that is further destabilized by e , 
however with the oimhd effect included the KBM appears only at higher (3 e (outside the 
plotted j3 e range). 
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FIG. 6: Imaginary and real parts of eigenvalues (7, co r ) as functions of /3 e for c p = 1 (solid lines) 
and c p = (dashed black lines). NSTX (a and d), ASDEX-UG (b and e), and JET (c and f). The 
vertical lines represent the base parameters color coded similarly to their respective machine. 



As the collisionality is suggested as an important parameter in driving the MTMs unsta- 
ble, here we examine the dependence of the MTM mode characteristics on this parameter. 
In figures [5] (b,d) the linear growth rates and the real frequencies of the most unstable modes 
are shown as functions of the collision frequency u e i for the different tokamaks. For the ex- 
perimental values of v ei (marked with vertical lines) the most unstable mode is found to be 
an MTM in all the three considered tokamaks. In NSTX, the growth rate increases with 
colMonality a 8 expected si nce colMonality is one of the instability drives fl Q, however, 
a further increase in v ei above the experimental value (almost doubled) the MTM growth 
rate does not increase significantly. In the ASDEX-UG the growth rate shows a decline as 
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collisionality increases. This trend is not surprising since previous studies have shown that 
the growth rate of the MTM has a non-monotonic dependence on the collisionality UEJ. As 
the collisionality further increases beyond the peak value, particles are very much scattered 
by collisions, and therefore preventing the formation of the current layer necessary for the 
MTMs to become unstable. Hence, as seen in figure [5] (b) we expect that this is the case 
for the ASDEX-UG case. For JET we also observe a (gentle) non-monotonic trend, and 
our base value of the collisionality seems to be positioned near the peak value. In all three 
machines no unstable MTMs were found for the collisionless case, i.e. u ei = 0, and under 
this condition the se p lasmas are found to be TEM unstable. Using the eigenvalue solver 
method in gyro 16] we followed the root corresponding to the MTM instability towards 
smaller u e i, and the mode is completely stabilized in the collisionless limit. The values of 
v ei jbj r corresponding to the maximum growth rate are rather different between the spherical 
and the conventional tokamak cases; in NSTX it is 2.05, while it is 0.55 and 0.51 in the 
AUG and the JET cases, respectively. 

In the literature the electron temperature gradient is suggested as one of the instability 
drives for MTMs js, 9]. This has been confirmed in previous numerical studies 4-6]. Here we 
compare the role of electron temperature gradient in destabilization of the MTMs between 
the three considered machines by performing a scan over the a/L^ e parameter. The results 
of this scan are shown in figures [7] (a,c) where the linear growth rates and the real frequencies 
of the most unstable modes are illustrated as functions of a/L^e- As seen in this figure, a 
finite value of a/Lx e is necessary for the destabilization of the MTMs in all three machines, 
however a clear difference is observed in the variation of MTM growth rates with a/L^e 
between spherical and conventional tokamaks. For the NSTX case, there is a clear increase 
of the MTM growth rate with an increase in a/Ly e , while for ASDEX-UG and JET the 
MTM growth rates show weaker and a non-monotonic dependence on ajL^ e . By further 
increase in a/Lxe the most unstable mode switches from an MTM to ITG/TEM modes, 
corresponding to the last points in the ASDEX-UG and JET curves in Figs. [7J (a and c). 

The a/Lxe threshold for MTM instability is observed to be well below our baseline values 
(indicated by vertical lines), and it is lower for JET and ASDEX-UG than for NSTX. The 
growth rate of MTMs is found to be less sensitive to electron temperature gradient for 
both of the conventional tokamaks than in the spherical tokamak studied here. In order 
to investigate the reason for this difference we performed a similar scan over a/Lxe and 
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FIG. 7: Imaginary and real parts of eigenvalues (7, cj r ) as functions of a/L^ e (a and c), and a/L ne 
(bandd). Red solid lines: NSTX, blue dashed: ASDEX-UG, green dash-dotted: JET. The vertical 
lines represent the base parameters color coded similarly to their respective machine. 

set c p = to eliminate the olmhd stabilization effect. The results are shown in figure El 
Without ocmhd stabilization the MTM still remains the most unstable mode for the baseline 
parameters, and as a/Lxe increases no significant change is observed for NSTX and ASDEX- 
UG. However, for JET case three regions in the a/Lxe space can be distinguished. For low 
values of a/L^e the MTM is the dominant instability; its growth rate increases linearly with 
a/Lxe- At a/Lxe ~ 2.5 an ITG mode, which was previously stabilized by olmhd effects - 
see figures [7] (a,c), - becomes the dominant instability, which gradually transits to a TEM 
as a/LTe becomes very large. Therefore, from our observations the stabilizing effect of the 
oimhd parameter on the MTM mode is not significant, but it has an impact on the stability 
of ITG/TEM/KBM modes. The strong suppression of these modes therefore, allows for 
MTMs to remain the dominant instability for a wider range of a/ ' L^ e . The weaker and non- 
monotonic dependence of the MTM on the electron temperature gradient in ASDEX-UG 
and JET thus, can not be explained by this effect. 

One of the main differences between the plasma parameters in the three plasmas consid- 
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FIG. 8: Imaginary and real parts of eigenvalues (7, u) r ) as functions of a/Lxe with c p = 1 and 
c p = 0. Red solid lines: NSTX, blue dashed: ASDEX-UG, green dash-dotted: JET. 

ered appears in the a/L ne values; a strongly negative electron density gradient, correspond- 
ing to a hollow electron density profile, is observed in NSTX, while the density profiles were 
slightly peaked in the ASDEX-UG and JET plasmas. Thus, we have investigated the depen- 
dence of MTM linear growth rates on the a/L ne parameter. Quasi-neutrality is enforced by 
slightly varying a/L ni while keeping the impurity density gradients fixed to the base parame- 
ters given in table [TT1 Figures [7] (b-d) show the linear growth rates and the real frequencies of 
the most unstable modes versus a/L ne for the different machines. In all three machines, the 
MTMs linear growth rates exhibit a non-monotonic dependence on a/L ne parameter with 
maxima corresponding to slightly hollow electron density profiles (a/L ne ~ —0.5). Clearly 
a/L ne is not a strong and necessary drive for the MTMs, as there are finite MTM growth 
rates in all machines at a/L ne = 0. 

In all the machines the MTM is the most unstable mode over the range of —1.5 < 
a/L ne < 0.5, with a rather weak dependence on this parameter. For sufficiently high electron 
density gradient the dominant linear mode switches from a MTM to an ITG mode in all 
plasmas. Again we examined the impact of olmhd by comparing the results to corresponding 
simulations with c p = (black dashed lines in figure [9]). Similar to the results of the a/Lxe 
scans shown in figure El the effect of the olmhd parameter on the MTM mode itself is 
negligible, but it strongly stabilizes the ITG mode. As shown in figure [HI without qcmhd 
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effects in the cases of NSTX and JET the a/L ne threshold of the dominant ITG mode is 
reduced, but in ASDEX-UG case no significant change is observed. 



NSTX AUG JET 




-9.5 -1 -0.5 0.5 1 -9.5 -1 -0.5 0.5 1 -9.5 -1 -0.5 0.5 1 




FIG. 9: Imaginary and real parts of eigenvalues (7, cj r ) as functions of a/L ne with c p = 1 and 
c p = 0. Red solid lines: NSTX, blue dashed: ASDEX-UG, green dash-dotted: JET. 



IV. CONCLUSIONS 



We have investigated the onset and parametric dependences of the MTM instability in the 
core (ptor = 0.6) of a spherical (NSTX), and two conventional tokamaks (ASDEX-UG, JET). 
The quasilinear transport is computed using the gyrokinetic code GYRO in the flux-tube 
(local) limit. In confirmation with previous studies, we found that for the experimentally 
relevant plasma parameters the MTMs are linearly the dominant instability in NSTX and 
ASDEX-UG. Under typical JET baseline parameters considered here the MTMs are also 
found as the dominant linear instability. In NSTX and JET the maximum of the MTM 
linear growth rate is located at higher mode numbers than that for ASDEX-UG. Therefore, 
the previously discussed idea that the higher mode number MTMs are the characteristics 
of the spherical tokamaks, while lower mode number MTMs are the characteristics of the 
conventional tokamaks is not supported by our results. 

Parametric scaling of the MTM instability in the core of these plasmas, revealed that 
a finite level of u ei and /3 e are needed in order for MTMs to become unstable, and when 
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unstable, they can remain the dominant instability over a wide range in v ei and (3 e . The 
linear MTM growth rate seems to be only weakly dependent on /3 e , and exhibits a weak but 
non-monotonic dependence on collisionality, v e i. By neglecting the collisionality or j3 e effects, 
the ITG/TEM modes appear as the most unstable modes in all three studied machines. 

A strong dependence for the growth rate of MTMs on electron temperature gradient is 
found in NSTX, while for ASDEX-UG and JET the MTM growth rate is found to be less 
sensitive to this parameter. The MTM growth rate significantly increases as a/L^ e increases 
in NSTX, but in ASDEX-UG and JET a weak and non-monotonic dependence on a/Lj<e 
is found. These results indicate that while a/LT e is a fundamental drive for the MTMs in 
these plasmas, it can contribute to the stabilization of the mode as well; this non-trivial 
behavior is more pronounced in the ASDEX-UG and JET plasmas. 

Similar trends are observed in all three machines when scanning for the electron density 
gradient, where the MTMs linear growth rates show again a non-monotonic dependences on 
a/L ne parameter with peaks located in the negative a/L ne region corresponding to slightly 
hollow electron density profiles. 

We have investigated the impact of a finite ocmhd on the onset and characteristics of 
MTM instability in various parametric scans, and we have observed that the stabilization 
of oimhd parameter on the MTM mode itself is not significant in all machines. However, its 
impact on the ITG/TEM unstable modes can result in a strong suppression of these modes 
allowing for the MTM to remain the dominant instability for a wider range in considered 
parameters, i.e. (3 e , a/Lxe and a/L ne . The non-monotonic and weaker dependence of the 
MTM on the electron temperature/density gradient however, can not be explained by cumhd 
effects. 

In the studied cases the MTM drives mostly electron heat transport through 6 An fluctu- 
ations, while other transport channels and contributions from 5<fr and 8B\\ are significantly 
smaller. For the ion heat flux and particle fluxes, the transport driven by MTMs are neg- 
ligible, however, if ballooning modes like ITG/TEM/KBM are also present, even as sub- 
dominant modes, these fluxes are mainly driven by these modes and therefore the overall 
ion heat and particle transport may not be negligible. However, these findings are based on 
linear analysis and further non-linear studies are needed for their confirmation. 

Finally, we would like to stress that the cases which we have considered have been chosen 
because they were unstable to MTMs. However, particularly large aspect ratio tokamaks, 
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being unstable to MTMs should not be considered a generic property. 
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